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ABSTRACT 



We show that dynamical relaxation in the aftermath of a galactic merger and the ensuing formation 
and decay of a binary massive black hole (MBH), are dominated by massive perturbers (MPs) such 
as giant molecular clouds or clusters. MPs accelerate relaxation by orders of magnitude relative to 
2-body stellar relaxation alone, and efficiently scatter stars into the binary MBH's orbit. The 3- 
body star-binary MBH interactions shrink the binary MBH to the point where energy losses from 
the emission of gravitational waves (GW) lead to rapid coalescence. We model this process based 
on observed and simulated MP distributions and take into account the decreased efficiency of the 
star-binary MBH interaction due to acceleration in the galactic potential. We show that mergers of 
gas-rich galactic nuclei lead to binary MBH coalescence well within the Hubble time. Moreover, lower- 
mass binary MBHs (< 10® M©) require only a few percent of the typical gas mass in a post-merger 
nucleus to coalesce in a Hubble time. The fate of a binary MBH in a gas poor galactic merger is less 
certain, although massive stellar structures (e.g. clusters, stellar rings) could likewise lead to efficient 
coalescence. These coalescence events are observable by their strong GW emission. MPs thus increase 
the cosmic rate of such GW events, lead to a higher mass deficit in the merged galactic core and 
suppress the formation of triple MBH systems and the resulting ejection of MBHs into intergalactic 
space. 

Subject headings: black hole physics — clusters — galaxies: nuclei — stars: kinematics — giant 
molecular clouds 



1. INTRODUCTION 

There is compelling evidence that massive black 
hole s (MB Hs) ex ist in the c e nters of most galax- 
ies (iFerrar esc & MerrittI 120001 : iGebhardt efal] 120031 
[Shields e"t~al.i i2003f l. It is beheved that galaxies grow 
by successive mergers, during which the two MBHs sink 
to the center of the newly formed galaxy by dynami- 
cal friction and form a "hard" binary MBH (BMBH) 
(jBegelman et al.lll980D with a semi-major axis of 



ah = [Q/(l + Q)']n.(Mi2)/4, 



(1) 



where M12 — Mi + M2 is the mass of the binary, Q = 
M2/M1 < 1 is the mass ratio and rfi{Mi2) is the radius 
of dynamical influence of the BMBH ^, where typically, 
ah ^ 1 — 10 pc. After the BMBH hardens, it continues 
to shrink by losing energy and angular momentum to 
stars, which are ejected fr om the system (the "slingshot 
effect" , ISasla w et aLlllQT^ or to gas with which it inter- 
acts dynamically. Once the separation further decreases 
by 2-3 orders of magnitude, the BMBH rapidly decays 
by the emission of gravitational waves (GWs) until the 
two MBHs coalesce. 

In mergers induced by interactions with stars, the 
merger timescale depends on the rate at which new 
stars are supplied to BMBH-crossing orbits ( "loss-cone" 
trajectories). Simulations show that this supply rate 



is typically not high enough; the BMBH stalls be- 
fore reaching a small enough separation for efficient de- 
cay by GW emissio n, and fails to coal esce in a Hub- 
ble time, tn (e.g. Berczik et al.l l2005l see review by 
iMerritt fc Milosavlie vic 2005^^ This result appears to 
contradict the circumstantial eviden ce that most galac- 
tic nuclei contain only a single MBH (jBerczik et al.ll2006l : 
IMerritt fc Milosavlievid [20051 ) . and furthermore, it im- 
plies few such very strong GW sources, which future GW 
detectors, such as the Laser Interferometric Space An- 
tenna (LISA), hope to detect. 

Several mechanisms were sugg ested as means of ac cel- 
erating BMBH coalescence (see IMerritt et al.l l2007l for 
a recent overview and discussion), either involving in- 
teractions with stars ("dry mergers") or with gas ("wet 
mergers"). These include re-ejection of stars that had 
a previous interaction with th e BMBH but were not 
ejected out of the galactic core ([Milosavlievic fc Merrittj 
120031 iBerczik e t ah '2005*): BMBHs embedded in dense 
gas (e.g. iIva noy__ct_al. ,.1999; ,Esc_aIaL_et_al,, .200j; 
iDotti et al.l I2007D: interactions of the BMBH with 
a third MBH (iMakino fc Ebisuzakil 11994 iBlaes et all 
I2002t llwasawa et al.l 120061 ): accelerated BMBH coales- 
cence due to accelerated loss-conc replenishment in non- 
axisymmetric potent ial ( Y u 2002; Berczik et al. 200(1) 
or in a steep cusp (jZieri [20061 I2007D . It is stiU un- 



clear whether these mechanisms are efficient enough, 
or whether they occur commonly enough to solve the 
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rcqmres the BMBH to be embedded m a very dense 



The William Z. & Eda Bess Novick career development chair 
^ D efined here as the ra dius that encloses a stellar mass of 
2Mi2 IIMerritt fc Szelll |2006| ) . The threshold semi-major axis for 
a hard BMBH is sometimes defined in terms of of cr^, the typi- 
cal velocity dispersion in the center, a/^ = GMiA'/2/i/4(T^, where 
fi = M1M2/M12 is the reduced mass. However, this is ill-defined 
since cr^ usually varies with distance from the BMBH. 



massive central gas concentration. However, it is un- 
known whether such amounts of gas exist there. For 
example, the central ~ 2 pc of the Galactic Center 
(GC) are gas-depleted (jChristopher et all l2005f l. and 
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othe r galaxies also show ce ntral gas cavities in their nu- 
clei (|Sakamoto et alJ [19991 ). The gas may also be dis- 
persed by the accreting BMBH before the merger is com- 
pleted (jMerritt fc Milosavlieyid I2005D . and ma y not be 
efficient for mino r mergers ("Esc ala et al.l 1200 5: but sec 
iDotti et al.ll2007l ). It is likewise unknown whether the 
non-axisym metric potential assu med by the dry merger 
scenario of iBerczik et alj ()2006[ ) is generally present in 
the post-merger galaxy on the relevant scales. Even if 
that is the case, actual demonstration of rapid BMBH 
coalescence still awa its future A^-body s imulations with 
reahstically high N (jBerczik et al.ll2006| ). 

Here we explore another possibility, which is likely to 
apply in most mergers: BMBH coalescence driven by 
massive perturbers (MPs), such as giant molecular clouds 
(GMCs) or stellar clusters, in the post-merger galaxy. 
MPs accelerate relaxation and scatter stars into the 
BMBH orbit at high ra tes. Efficient relaxation by MPs 
was first suggested by iSpitzer fc Schwarzschildl ()1951l . 
|1953[) to explain stellar velocities in the galactic disk. 
MPs remain an important component in modern models 
of galactic disk heating (see e. g. Villumscn 1983, 1985; 
Lacevlligsl : iJenkins fc Binnevlll990; .Hanninen fc Flvnnl 
20021 and references therein). A similar mechanism was 
proposed to explain the spatial diffusion of stars in the in- 
ner Galactic bulge (Kim fc Morris 2001). In addition to 
dynamical heating, efhcient relaxation by MPs was sug- 
gested as a mechanism for loss cone replenishment and 
relaxation, both in the context of scattering of Oort cloud 
comets to the Sun ()Hillslll981l : [Bailevlll983l ) a nd the scat- 
terin g o f stars to a MBH in a galactic nucleus (|Zhao et al.l 
I2002D . IZhao et al.l (120021 ) also noted the possibility of 
increased tidal disruption flares and accel erated MBH 
binary coalescence due to MPs. Recently, iPerets et al.l 
(|2007l ) (Paper I) studied in detail MP-d riven interactions 
of single and binary stars with a single MBH. 

In this study we apply the methods developed in Pa- 
per I to investigate MP-driven interactions of stars with 
a BMBH, and the consequences for BMBH coalescence. 
We explore different MP populations and merger sce- 
narios based on the available observations and simula- 
tions, and estimate the BMBH coalescence rate. We also 
discuss additional implications: the mass deficit in the 
galactic cores and the suppression of triple MBH forma- 
tion in galactic mergers. 

This paper is organized as follows. We begin with a 
general overview of our calculations and the new results 
(Sj2]). In Sj3] we briefly summarize the physics of MP- 
driven loss-cone replenishment, which are derived in de- 
tail in Paper I. The observations and theoretical predic- 
tions of MPs in the inner hundreds pc of galactic nuclei 
are reviewed in SQ] and used to construct the MP models 
used in our calculations (fj5]) . In ^l6]we briefly review the 
dynamics of BMBH mergers; a detailed technical discus- 
sion is presented in appendices lAl andlBl We then present 
our procedure for modeling the dynamical evolution of 
the BMBH coalescence under various assumptions in Sj6] 
and analyze the results of our calculations in fjT] We ex- 
plore their implications in fJ5]and discuss and summarize 
our results in ^ 

2. OVERVIEW OF THE CALCULATIONS AND NEW 
RESULTS 



This paper presents a first detailed study of the im- 
pact of MPs on BMBH mergers. In this section we give 
an overview of the methods, calculations and modeling 
that we use to study the dynamics of BMBH mergers by 
MP-induced slingshot events, and to estimate the merger 
timescale and the associated mass ejection from the nu- 
cleus. These are described in detail in ^31-^ and in the 
appendices. We then briefly list the new results derived 
here, and their implications. These are described and 
discussed in detail in SjTj-SjSl 

MPs shorten the BMBH merger timescale by acceler- 
ating the supply rate of stars to the loss-cone (see 
The application of loss-cone theory to the BMBH merger 
problem requires a model of the MP properties: their 
mass, size, number and spatial distributions. We review 
the available observations and numerical simulations of 
MPs in merging galaxies of different morphological types 
(21 tables [T] and [2]). We then construct a few generic 
models for various types of merger scenarios, with and 
without MPs, based on these data and on theoretical 
resuhs (gSl table O. 

To evolve the BMBH in time in the context of a given 
merger scenario (deflned by a BMBH mass ratio and 
a MP model, see we execute the following steps. 
(1) We establish the BMBH's initial semi-major axis at 
the stalling point, where 2-body stellar relaxation can 
no longer efficiently resupply the loss-cone. The stalling 
radius is estimated semi-analytically by the formula de- 
rived in appe ndix El which are based on the A'^-body 
simulations of iMerritH ()2006l ). (2) At each time step, 
we calculate the rate at which stars are scattered into 
the loss-cone by MPs from all radii taking into account 
the BMBH's instantaneous separation (Paper I; [2]). (3) 
We calculate the hardness of the encounter (the ratio 
between the encounter velocity and the BMBH orbital 
velocity) due to the acceleration of the infalling star 
in the galactic potential, as derived in appendix IB. 21 
This effect, which was neglected in previous studies, sig- 
nificantly decreases the energy extraction efficiency of 
the encounter. (4) We estimate the average amount of 
BMBH orbital energy extracted by the interaction with 
the stars scattered into the loss-cone. This is based 
on numerical e xperiments of isolated 3-body encounters 
(|Quinlanlll996l hereafter Q96), which are adapted in ap- 
pendix |B?T] to take into account the changing hardness of 
the encounter in the course of the BMBH evolution. (5) 
We include the energy extraction by GW emission, which 
increasingly dominates the last stages of the merger. (6) 
We update the BMBH semi-major axis according to the 
extracted energy by both the slingshot effect and GW 
emission. (7) We end the dynamical calculation when the 
emitted GW power exceeds that extracted by dynamical 
interactions, and the BMBH enters the final rapid GW- 
induced decay phase leading to coalescence. 

Our new results and their implications are as follows. 
(1) The inferred properties and numbers of MPs in post- 
merger galactic nuclei are high enough to effectively re- 
duce the relaxation times there by orders of magnitude. 
Consequently, BMBH decay due to MP-induced dynam- 
ical interactions with stars is efficient enough in most 
mergers to lead to coalescence within a Hubble time (typ- 
ically within 10®~^ yrs), thereby solving the "last par- 
sec" problem. The rates of GW emission from BMBH 
coalescence should thus trace the galactic merger rates. 
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(2) Since the MP-induced dynamical phase of BMBH 
mergers is rapid, we predict that most BMBHs will be 
observed either electromagnetically at separations larger 
than the hardening radius (> Ipc), or by GW emission 
at much smaller separations (typically <C 1 pc), during 
the GW-dominated decay phase. (3) In most cases the 
BMBH coalescence takes less than the typical time be- 
tween major galactic mergers. Therefore the formation 
of triple MBHs through multiple galactic mergers is un- 
likely, and ejection of MBHs from unstable triples should 
be rare. (4) The ejection of stars by the slingshot effect 
removes mass from the galactic core ("mass deficit"). 
The more the BMBH can decay, the larger the total 
ejected mass. The spatial scale of the deficit reflects the 
origin of the scattered stars. We point out that the mag- 
nitude and spatial extent of the deficit can be used to 
discriminate between merger mechanisms. Specifically, 
the MP-induced BMBH merger mechanism predicts a 
larger mass deficit over a larger spatial extent than that 
due to stalled mergers by stellar relaxation alone. 

3. LOSS-CONE REFILLING BY MASSIVE PERTURBERS 

In Paper I we present a detailed quantitative analysis of 
the MP-induced resupply of stars to nearly radial , "loss- 
cone" orbits, which bring them within some threshold 
distance q from the central mass (MBH or BMBH) , where 
they undergo a strong interaction with it ( "event" ) and 
are thereby removed (scattered or destroyed). We show 
that the resupply rate by MPs is orders of magnitude 
faster than that by stellar 2-body relaxation alone. This 
translates to an accelerated rate of close interactions with 
the central object. Here we present a brief qualitative 
summary of the results of paper I and their implications. 

A star with orbital energy in the range {E, E+dE) on 
a nearly-radial loss-cone orbit with angular momentum 
J < Jic will reach the MBH and be removed in about 
a single dynamical (or orbital) time P{E). When the 
resupply rate of such stars is slower than the rate at 
which they are removed {^1/ P), the phase-space region 
(dE, J < Jic) is nearly empty of stars. Conversely, when 
the resupply rate is higher than 1/P, that phase-space 
region is nearly full, and its phase-space distribution is 
nearly isotropic. At that point the rate at which stars 
with energy E interact with the central mass reaches its 
maximal value; further scattering can not increase the 
rate. When the resupply of stars is driven by stellar 
2-body relaxation, tightly bound regions of phase space 
(high E, small typical r) are empty, since the angular size 
of the loss-cone is large, the period is short and removal 
is fast. Conversely, loosely bound regions (low E, large 
typical r) are full. Most of the contribution to the total 
resupply rate comes from stars with energies near the 
transition between the empty and f ull loss-cone regimes , 
where P ^tr, the relaxat ion time (|Frank fc ReedlTOTd 
iLightman fc Shapirolll977[ ). This simple picture changes 
when MPs dominate relaxation, because their distribu- 
tion does not necessarily follows that of the stars, and 
their properties (mass, size) may change with distance 
from the center. 

MPs of mass Mp and space density Up dominate dy- 
namical relaxation over scattering by stars of mass M^, 
and space density n^, when the ratio of the 2nd moments 
of the mass distributions satisfies //2 = ripMp j n^^Ml > 
t*/t^^ ^ 1. The quantity /i2 thus expresses the MP 



enhancement factor of the relaxation timescale, up to 
the ratio of the Coulomb logarithms for star-star and 
star-MP scatterings, which takes into account the ex- 
tended size of the MPs . This can be shown by con- 
sidering first close encounters at the "capture radius" 
Tc ~ GMp/v^, where v is the typical relative velocity. 
The "nuS" rate of such encounters with a test star is 
then ~ nvrl (x UpMp/v^. Integration over all en- 
counter distances further decreases the relaxation time 
by a Coulomb logarithm factor that depends on the 
size of the system and, if Rp > Tc, also on the 
size of the MP"^. The impact of fast, MP-induced re- 
laxation on the rate at which stars enter the loss-cone 
can be incorporated into standard loss- cone theory (e.g. 
ILightman fc Shapirolll977l : lYoungjflQTTl ) by replacing the 
relaxation time due to scattering by stars with that due 
to MPs, with the only modification being the separate 
treatment of rare scattering events. 

The effect of MPs on the BMBH loss-cone depends on 
the typical distance from the center and on the BMBH's 
stage of dynamical evolution. Far from the BMBH, MPs 
fill the loss-cone and drive the resupply rate to its maxi- 
mal value, while closer to the BMBH, even MPs can not 
fill the loss-cone, but they still increase the supply rate 
by a large factor ^ ii2- Closer still, the tidal field of the 
BMBH gradually limits the size and mass of the MPs 
to the point where stellar relaxation re-emerges as the 
dominant relaxation mechanism (see Fig. 3 in Paper I). 
The demarcation between the full and empty loss-cone 
regions changes with time as the BMBH decays, the size 
of the loss-cone decreases and the empty loss-cone region 
contracts. Irrespective of these details, the overall effect 
of MPs is to dramatically increase the loss-cone resup- 
ply rate over what would be expected from slow stellar 
2-body relaxation alone. 

We calculate the loss-cone resupply rate with some sim- 
plifying approximations. We assume a spherically sym- 
metric distribution, the ansatz E <-> GM{< r)/2r and 
the Keplerian expression for the orbital period P{E). We 
also make the conservative assumption that one class of 
massive perturbers dominates relaxation, and neglect the 
contributions from other mass components (e.g. stars). 
The full expressions for the supply rates of stars deflected 
from a typical radius r in the empty loss and full loss 
cone regimes (§2 in paper I, ) can be respectively ap- 
proximated by 



dPe N,{<r) ATf q N4<r) 

dlogr \og{r/q)tr ' dlogr r P{r) 



(2) 



where Ni,{<r) is the number of stars enclosed within r, 
and where the dependence on the MPs properties enters 
through the relaxation time, tr — oc /i2i*. The total 
loss-cone refilling rate is then given by integrating over 
the contributions from all radii, taking into account the 
fact that the refilling rate does not exceed the full-loss 
cone rate. 



dTejq) dTfjq) 



dr 



dr 



dr . 



(3) 



^ For example, on a scale of r~50 pc where the velocity disper- 
sion is (T~ 0(100 kms~^) and for MPs of typical size Rp ^ 5 pc, 
At-^ro-^/GA/*, Aup^r/Rp and log A*/ log Amp ~ 8. 
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In the context of BMBH decay, the rate (Eq. [3]) implicitly 
depends on time through the evolution of q = a(t), and 
the rate, in turn, affects the evolution of a through the 
BMBH evolution equation (Eq. [7] below). The BMBH is 
evolved in time by iterative numerical integration of the 
coupled rate and evolution equations. 

The effects of MPs are significant in situations where 
perturbations by stars alone are not efficient enough to 
refill the loss-cone on the orbital timescale. MPs are most 
effective for large-q processes, such as close interactions 
between binaries and a single MBH (where q is the tidal 
disruption radius of the binary), or interactions of single 
stars with a BMBH (where q is the BMBH's semi-major 
axis). This is because the larger q, the lower the empty 
loss-cone refilling rate relative to that needed to com- 
pletely fill the loss-cone (i.e. Pe/Pycx (r/g)/log(r/q) is a 
decreasing function of q for r/q^l, Eq. 

4. MASSIVE PERTURBERS IN GALACTIC NUCLEI 

The space density of MPs is much smaller than that 
of stars, so to dominate relaxation (/i2 3> f ) they must 
be significantly more massive. Here we consider only 
MPs with masses Mp > W^Mq, such as stellar clusters 
and GMCs gas clumps. Intermediate mass black holes 
(IMBHs) could be very effective MPs, but these are not 
considered here since it is still unclear whether they ac- 
tually exist. A summary of the observed properties of 
MPs and those derived from simulations is presented in 
tables [T] and H 

The lifespan of single MPs is affected by many factors, 
such as collisions with other MPs; the central galactic 
tidal field, which can disrupt them if they approach the 
center on an eccentric orbit or sink there by dynamical 
friction; GMCs can also be dispersed by strong radiation 
from stars and by supernovae. Thus, individual MPs are 
not expected to survive much longer than a few local 
dynamical times, a time substantially shorter than the 
galactic merger timescale (roughly the galactic dynami- 
cal time), which in turn may be much shorter than the 
time required for BMBH merger. However, MPs do not 
have to exist individually much longer than a local dy- 
namical time to have a strong effect on relaxation. It 
is sufficient that the overall MP population maintains 
a steady state, by the continuous formation or supply 
of new MPs to replace those destroyed. There is both 
observational and theoretical support for the long-term 
persistence of large populations of MPs in galactic nuclei. 

Observation of the GC indicate the steady state ex- 
istence of a large population of short-lived MPs over a 
Hubble time. As reviewed in detail in Paper I, observa- 
tions of dense gas in the central ~ 100 pc of the GC reveal 
Np ~ 100 GMCs with a t ypical mass of (Mp) ~ 10^ Mq 
(table[T]) (|Oka et al.l2001h . GMCs and young stellar clus- 
ters are stages in the path of star formation, and so the 
star formation rat e can b e used to estimate the MP mass 
supply rate. Figer et al] ()2004f ) show that the star for- 
mation history in the central projected 30 pc of the GC is 
well described by continuous star formation over 10 Gyr 
at a rate of 0.02 Mq yr^^. Extrapolated out to 100 pc in 
the n^, ^ stellar distribution of the inner bulge, this 
corresponds to dM^,/dt^0.05 Mqjt^^. Since the mean 
star formation efficiency (fraction of mass turned into 
stars) is on average very low, f^, ^ few x 0.01 (jMvers et al.l 
Il986f) . the star formation rate translates to an MP mass 



supply rate of dM/dt - {dAh/dt)/f^ - 0{lMQyr-^) 
and MPs formation or supply rate of 



dM/dt 



0.01 / V 105 M, 



'0 



(4) 



The mean lifetime of such GMCs is then 

tp^^^^IxW yi ^^^^j [io^MqJ [lOO 

which is comparable to the dynamical time scale at ~ 100 
pc. We therefore conclude that the observed MPs in the 
GC, together with the inferred star formation rate and 
history, are fully consistent with such a steady state MP 
population persisting over a Hubble time. 

Simulations of gas in galactic nuclei also show highly 
inhomogeneous quasi-steady state conditions, with a very 
broad mass spectrum extending over ^ 7 or ders of magni- 
tude (|Wadall200lUWada fc Nonna^[200lf) . The nuclear 
MP population in the Galaxy is probably representative 
of that in nuclei of other spiral galaxies, where similar 
amounts of dense ga s are observed (Kodae t al. 200^; 
ISakamoto erallll999l : lYoung fc ScoviilillQQlh . Further- 
more, estimates of the gas supply timescales for star 
bursts in other galaxie s also suggest c ontinuous gas sup- 
ply over 10*^ - 10^ yr (|Leithereii l200lV 

The MP contents in the nuclei of post-merger 
galaxies are expected to be yet larger than in the 
Galaxy, which does not appear to have undergone 
a recent major merger. This is indeed observed in 
UL IRGS jS anders fc Mirabel' n^96^, in merging galax- 
ies (lEvans et al. 2002; Sakamoto ct al. 2006; CuUcn et all 



2007), and also seen in simulations ([Barnes fc Hernquis 



19921) 



We now turn to a detailed discussion of the observed 
and inferred properties of MPs in galactic nuclei of differ- 
ent types. These form the basis of the MP models used 
in our numerical analysis of BMBH mergers. 

4.1. Massive perturhers in spiral galaxies 

The ~ 100 GMCs observed in the GC, with masses in 
the range IO-^-IO^Mq, contain few x 0.01 of the total 
dynamical mass on the few x 100 pc scale and a few x 0.1 
in the central ~ 100 pc (see Paper I for an extended dis- 
cussion of the properties of MPs in our GC). In contrast, 
the central ~2 pc of the GC contain negligible amounts 
of gas. Single molecular clouds cannot be resolved in 
the nuclei of other spiral galaxies, but the total fraction 
of gas and its distribution are usually quite similar to 
those observed in t he GC (e.g. Sakamoto ct al. 199^; 
ISawada et all 120041 : see review bv iHenkel et al.]ll99l[ l. 
CO observations show that the gas contains very dense 
large clumps that account for up to < 50% of the to- 
tal gas contents in these regions (jPownes et al.l 119931 : 
.Downes fc Solomon 1998,). 

In addition to GMCs, m any globular clusters (jFriell 
IT995I: lAshman fc Z"e^[l998h and open clusters may in- 
spiral into, or form in t he galactic nucleus in the course 
of their evolution (e.g. iGnedin et al.|[l999( ). For exam- 
ple, the Galaxy contains hundreds of ~ 10^ Mq open clus- 
ters and fewxlO^ MfTi globular clusters (jMeylan fc Mavoil 
fl99lt lMeilfT995h. Many more are observed in other 
galaxies ( Ashman fc Zepj|1998f ) . If some of these clusters 
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contain IMBHs, they will contribute to the MPs popula- 
tion eve n after the disruption of the host cluster is dis - 
rupted ()Ebisuzaki et all 1200 ll : iMiller fc Hamiltonll2002D . 
and will sink all the way to the center. 

4.2. Massive perturbers in elliptical galaxies 

The gas fraction in elliptical galaxies is t ypically 10 — 
100 times sm aller than in spiral galaxies dRupenl [l997l : 
lKnappiri999f ). However, in some elliptical galaxies it is 
comparable or even larger than that in spirals. Such gas- 
rich ellipticals are thought to have been formed recently 
in a m erger of two late type galaxies (e.g. IWiklind et al.l 
Il997f) . In particular, ultra- luminous infrared galaxies 
(ULIRGs, see review by Sanders & Mirabel 1996) have 
extreme amounts of gas, 10 — 100 times more than in 
the Galaxy, and can have as much or more mass in gas 
compared to the mass in stars. Elliptical galaxies may 
well be evolved merger products, where m ost of the dense 
gas in the core formed stars (e.g. Bender & Sagiia 1999). 
In that case, it is plausible that the main type of MPs 
would be the stellar clusters that were born of the GMCs, 
rather than the GMCs themselves. Observations of stel- 
lar rings and disks in the cores of elliptical galaxies indeed 
suggest that present-d ay stellar structures refle ct ear- 
her gaseous structures (iDownes fc Solomonl[l99l . This 
is also consistent with the fact that ellipticals have 
larger numbers of globular clusters than spirals, and 
that mergers are associated with the for mation of mas- 
sive clusters jAs hman & Zepf 1998; Zh ang fc FalU 119991: 
iKravtsov feG^dinl 120051 : lLarsenll2006f ). 

4.3. Formation of massive perturbers in galactic 
mergers 

Simulations of mergers of gas rich spirals indicate 
that > 50% of the total gas mass in both galaxies is 
driven i nto the central few x 100 pc of t he newly formed 
galaxy (|Barnes fc Hernauistlll99lUl996D . where it prob- 
ably forms massive clumps. In mergers of two gas- 
poor ellipticals, stellar clusters may play a similar role. 
Many of the newly formed stell ar clusters will proba- 
bly su rvive in the merged nucleus (|Portegies Zwart et al.l 
I2OO2D . In addition, many old globular clusters will fall 
directly into th e nucleus in the course of the merger 



Gnedin fc Prieto 2 006). or sink in by dynamical friction 



Capuzzo-Dolcett al 1 1 9931 ) . While most will probably be 



disrupted (O. Gnedin, priv. comm.), a si gnificant frac- 
tion could survive (e.g. simulations by iMiocchi et al.l 
120061) . This central accumulation of young and old stel- 
lar cluster could significantly shorten the relaxation time. 
Further simulations are needed to address these questions 
quantitatively. 



5. MODELING MASSIVE PERTURBER-DRIVEN BMBH 
COALESCENCE 

Based on the observations and simulations described 
above, we formulate three representative merger scenar- 
ios that include MPs, and compare them to a merger sce- 
nario where only stellar 2-body relaxation plays a role. 
The model parameters are listed in table [3] 

The major merger scenario consists of a Q = 1 merger 
of two gas-rich galaxies. It is assumed that the merger 
triggers a large gas inflow to the center, increasing the 



amount of gas there to ~ 1/2 of the total dynamical mass 
(~ 5 times more than presently in the center of the Milky 
Way; the mass of the cold gas in post-merger galaxies can 
be even higher, but we take into account only the dens- 
est regions that correspond to the more massive MPs). 
It is further assumed that the MPs are similar to mas- 
sive GMCs in our GC, that they have a power-law mass 
function, dNp/dMp oc Mp/^ with (3 = 1.2 { see MP model 
GMCl in Paper I for details), and that their spatial dis- 
tribution is isotropic"*. 

The minor merger scenario consists of a Q = 0.05 
merger between a large, massive gas-rich galaxy and a 
much smaller galaxy, which only slightly perturbs the 
large galaxy and triggers only a moderate gas inflow to 
the center. It is assumed that the nuclear gas mass is 
~ 1/3 of the total dynamical mass (~ 1.5 times more 
than presently in the center of the Milky Way). The MP 
properties are the same as in the major merger scenario. 

In the elliptical merger scenario we attempt to model 
a Q = 1 merger of two equal mass gas-poor elliptical 
galaxies. We assume that the MPs are mostly stellar sys- 
tems such as clusters or spiral structures. Lacking secure 
observati ons, we model the MPs after results fr om sim- 
ulations (|Li et al.l [200I iPrieto fc GnedhH [200I) . These 
simulations show that both the total cluster birth-rate 
and the massi ve cluster bir th-rate peak at the center 
of the galaxy (|Li et al.ll26oi) . We assume that the MP 
mass fraction is 0.2 of the total dynamical mass and that 
the cluster mass function is a power-law with f3 — 2 
for 10^ < Mr, < W Mq^, following the results of 



IPrieto fc GnedhJ ^OM). 

Finally, we consider, for comparison, a model that as- 
sumes that relaxation in the post-merger galaxy is due 
to stellar 2-body interactions only. 

In our calculations we assume that the stellar distri- 
bution over the entire relevant distance range can be ap- 
proximated by a singular isothermal stellar distribution 



p(r) 



27rG'r2 ' 



(6) 



where the velocity dispersion CToo, and hence the nor- 
malizatio n, is determined by th e empirical M,/a rela- 
tion fe.g. iWang fc Merritti [200l . The MP distribution 
is assumed to follow the stars, down to a minimal radius 
'"MP, where the MPs are destroyed either by the cen- 
tral tidal field, the radiation of the accreting BMBH, or 
the outflows associated with the accretion or star forma- 
tion triggered by the merger. The exact value of tmp 
is uncertain, since the processes involved in the destruc- 
tion of the MPs are complex, and t he inner cusp may 
be fl attened due to previous mergers (jMilosavljevic et ahl 
l2002ai) . Here it is assumed that tmp = 2rh ^. In the 
GC, tmp ~ 4-8 pc. At r > tmp our assumed minimal 
CMC mass of 5 x 10"* M p) is consiste nt with observations 
(jBacker fc Srameklll999t lOka et al.ir2001h . This minimal 
radius is probably a conservative estimate, since tran- 
sient dense clumps and dense cluster cores can survive 
even at smaller distances, as seen in observations in our 

* While the geometry of central molecular zone of the Galaxy is 
flattened, its height of few X 10 pc implies that it is nearly isotropic 
of the scale of interest of ~ 100 pc. 

^ Note that the M./cr and M./M;, relations {Mi, is the mass of 
the bulge, with typical length scale rj) then imply that r^po^ri,. 
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Table 1 

Observed and simulated properties of massive perturbers 

MP type Mp (Mq) Mass Profile (Mp)^^^ (Mq) Rp (pc) References 

GMCs in tlie GC lO"* - 10** Power law (/3 = 1.2) 4 X lO'^ 5 lOka et al.l 1I2OOII ') : [Gusten fc Philippl 1I2OOI) : 

Perets et al. (20071 

Young clusters in the GC 10^ - 10^ Power law (/3 = 1.2) 3x10* 1 IFiger et al.l ||T999| . [2003) : IMaillard et al.l 1I2OOI) : 

Borissova ct al. (2005): Perets et al. (2007') 
Globular clusters in the Galaxy 10^-^ - 10^'^ Log normal 1.9 X 10^ 5 Mandushcv ct al^ (.199lh 

Young clusters in galaxies 10"'-^ - lO^'^ Power law (/3 = 2) 4.3x10^ 3 Zhang fc Fall ( 19991) : [Kravtsov fc GnedinI ||2005|') : 
Larsen (2006) 



Table 2 

Mass fraction of observed and predicted 

MASSIVE perturbers IN GALACTIC NUCLEI 

Galaxy MP type Afp°VMdyn References 

Milky Way GMCs few X 0.1 lOka et al.l jgOOj) ; [Husten fc Philip pl pOOl) 

Clusters 10"'' 'Figcr ct al. (1999, 2002); Maillard ct al. (200?); 'Borissova ct al.' ("200^) 

Spirals GMCs 0.1-0.3 iKoda et al., 1,2005, '): .Glisten fc Philipp 12004) : .Sakamoto et al.i (il99a ): lYoung fc Scovilld 1IT991I) 

EUiptical mergers (Sim.) Clusters 0.2" (Li et al. 2004 : Prieto fc Gnedin 2006) 

Ellipticals (Obs.) GMCs lO'Mo-^ IRupen) lfT997l) : IKnappI ||T999|) 

ULIRGs GMCs 0.3-0.6 ISandcrs fc Mirabel (19^) 

Merger (Obs.) GMCs 0.3-0.6 |Evans ct al. (2002); Sakam oto eTaTI 1I2OO6I ): ICuUen et al.l i200% 

Merger (Sim.) GMCs 0.3-0.6 [Barnes fc Hcrnquist (1992) 

" Assumed, based on simulations. See text. 



GC; (Christopher et al."'2005'), a nd found in simulati ons, 
fWada fc Norm an 2001; Porteg ies Zwart et aLl[200l ). 

6. BMBH MERGER DYNAMICS 

A BMBH m erger progresses through three stages (See 
lMerrittll2006[ ). (1) Gradual decay by dynamical friction 
to the point where the separation between the two MBHs 
is ri2 ~r/i(Mi). (2) Formation of a bound Keplerian pair, 
when ri2 <rh{Mi), through rapid decay, initially by dy- 
namical friction on M2 and later by the slingshot effect. 
This is followed by a slow-down of the decay when a^ah 
and stalling, unless the the loss-cone is replenished by 
a process more efficient than diffusion due to 2-body re- 
laxation. (3) Ultimately, the BMBH orbital decay rate is 
dominated by GW emission, leading to final coalescence. 
The operational definition of the stalling separation Os 
at time tg is the point where the decay rate sharply de- 
creases. Typically as'--^0(a?i) (see appendix R|l . 

The slingshot effect occurs when g, the periapse dis- 
tance of the star from the BMBH center of mass, is of 
the order of the BMBH semi-major axis a. Such stars 
are ejected and lost from the system, either directly or 
after several repeated interactions with the BMBH, and 
on average extract energy AE{q) from the BMBH. The 
evolution of the BMBH energy, or equivalently, the de- 
crease in a, is given by 



d /GM1M2 
di \ 2a 



'dr 

d^ 



AE{q)dq = r(a) {AE){a), 



(7) 

where T{a) is the supply rate of stars that approach 
the BMBH on orbits with q < a, and {AE) cx a"^ 
is the appropriately we i ghted mean extracted energy 
(|Milosavlievi6 fc MerritH l2003t iMerritt fc Milosavlievid 



see detailed discussion in appendix IB.ip . ft then 
follows from Eq. ^ that the dynamical decay rate in 
the two regimes scales as ddyn oc — a/ log(r/a) or oc— a^, 
respectively, so that in both cases the dynamical harden- 
ing rate decreases as a decreases. Note that in the hard 
BMBH limit (a — > 0), when the loss-cone is small and 
therefore full, d(l/a)/dt- const. (Q96). 

When the BMBH separation becomes small enough, 
the orbital decay rate due to GW emission, dcwi be- 
comes higher than the dynamical decay rate. We con- 
servatively assume circular BMBHs (eccentric BMBHs 
coalesce faster in the GW emission dominated phase). 
The de cay rate on a circular orbit due to the emission of 
GW is (lPeterslll964D 



64 G^^iMl2 



(8) 



which increases as a decreases. The time to decay to 
a = from an initial semi-major axis a is 



256 G3 /iMfj 



(9) 



Since ddyn decreases with a, while dew increases, there 
exists a transition BMBH separation, acwi such that 
ddyn(aGw) = aG'w(aG'w)- Once the BMBH shrinks to 
ogw, the coalescence is inevitable as long as tG"w(aGw) < 
Ih and as long as the BMBH remains unperturbed. The 
total time from the hardening semi-major axis ah to the 
coalescence is then 

tc = tdyn{ah flGw) + tGwi^GW " > 0) . (10) 

The dynamical decay timescale tdyn(ah ^ ^gw) is 
of the order of the time it takes the BMBH to inter- 
cept and interact with stars whose total mass equals 
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Table 3 

Massive perturber models in galactic mergers 

Merger model Q r/ry,'^ M*°VM*°*^ Mp{Mq) (3'' Rp (pc) 112" 

Major 1 2-30 1/2 5xl0"'-lxl0^ 1.2 5 3x10'' 

Minor 0.05 2-30 1/3 SxlO-'-lxlO^ 1.2 5 5x10^ 

Elliptical 1 2-30 1/5 1x10^-1x10^ 2 3 5x10^ 

Stars — 1-30 1 1 — 1 

" Spatial extent of the MP distribution, assuming Np{r) Qcr~^ . 

Assuming dNp/dMpOcMp^ 

" fi2 = Np (Mg) /Af^ (Af^) , where (M^) = J M'^{dN/dM)dM /N . 



its own , tdyn ~ Ml2/[M^,^{aGw)], where T is eval- 
uated at ttGW) where the rate is slowest^. This esti- 
mate neglects the possibility that a fraction of the stars 
are not ejected from the loss-cone, but return to inter- 
act again with the BMBH. This can further accelerate 
the decay, but is not enough i n itself to prevent stalling 
(|Milosavlievic fc Merrit"tll2003f) . 

In order to calculate the MP-induced coalescence time, 
it is necessary to compute both the rate at which stars 
are scattered by MPs into the BMBH, and the energy 
they extract from the BMBH . Based on results from 3- 
body scattering experini ents (jHillsl Il983t iQuinlanI 119961 : 
iSesana et al1[2fm l200l . we assume that a star whose 
periapse distance q from the BMBH's center of mass is 
smaller than the BMBH separation a, interacts strongly 
with the MBH and is then ejected out of the system. 
We omit the possibility of re-ejection, and we neglect 
soft scattering events (g > a) since these are inefficient 
in ex tracting energ y fro m the BMBH (see jSesana ct al. 
I2OO6I and appendix IB. ip . We thus obtain a conservative 
upper limit on the coalescence time. 

Beginning with a hard BMBH of separation a{t = Q) = 
Qs (appendix [Aj, we define the time-dependent loss-cone 
periapse as q — a{t) and calculate the loss cone rate ^{q), 
using the methods described in Paper I and summarized 
here in fj2l We follow the evolution of the BMBH sep- 
aration by numerically integrating the evolution equa- 
tion with small enough time steps suc h that da <C a 
(seepilosavljcvic & Mcrritt 2003 and Scs ana et ani2006l . 
I2OO7I for a similar approach) until the orbital decay is 
dominated by GW emission (Eq. [TU]) . which is effec- 
tively the coalescence ti me tc - In simplified notation, the 
evolution equation (Eq. IB15|) is 



dloga _ _2 



dt 



CM^ir, (11) 
dr 



where dV [dr is the differential loss cone replenishment 
rate and C is the mean value of the dimensionless energy, 
C= Myi^E 1 2Mi,Ei2, exchanged between the scattered 
star and the BMBH (see detailed derivation and numeric 
estimation in appendix lB.il C = 1 corresponds to the case 

" Every star that passes near the binary MBH extracts from 
it binding energy of order Mtei2, where £12 = Gii/2a is the 
specific energy of the BMBH, so that dE = -GhhM2/2a?da = 
{MirGfi/2a)r(a)dt. Integrating between ai^ 2> with r{a) ^ 

\{a/r)Nir{< r)/P{r)\^ (when the loss-cone is filled by MPs 

that orbit as close as r^p from the MBH) yields tdyn — 
Mi2/A/*r(aGw)- 




Fig. 1. — The dynamical decay times, t^yn of BMBHs from 
to oqw ^ function of the BMBH mass (the time to final GW- 
induced decay, fowi from a^w to a = is negligible compared 
to the initial dynamical decay phase). Different merger scenarios 
are shown (table [Sj: major mergers (solid line), minor mergers 
(dashed line) and elliptical mergers (dashed-dotted line). Without 
MPs (crosses, for Q = 1) the decay times are longer than the age 
of the universe (horizontal line) for all BMBH mass in this range. 



where the specific energy carried by the star equals twice 
that of the BMBH). The quantity C{a,r) depends on the 
hardness parameter of the encounter C = a{r)/Vi2{a), 
defined as the ratio between the typical initial velocity 
of the scattered star far from the BMBH, o'(r) and the 
orbital velocity of the BMBH, V12 = v^GMiz/a. An 
additional r-dependence is introduced by the accelera- 
tion of the star toward the BMBH by galactic poten- 
tial, which increases the relative velocity between the 
BMBH and the star at the point of encounter over what 
it would have been if the star fell toward an isolated 
BMBH (see appendix lB.2p and decreases the efficiency of 
the slingshot effect (Figure [3]). This non- negligible effect, 
taken into account here, was neglected in previous esti- 
mations of the BM BH coalescence times (jSesana et al.l 
[200l;[Quinia3[i99i). 

7. RESULTS 

Figure ^ shows the total decay time of BMBHs in 
the mass range M12 = lO^-lO^M© for different merger 
scenarios. Stellar 2-body relaxation cannot replenish 
the loss cone fast enough. In the absence of MPs, the 
merger proceeds in the empty loss-cone regime, where 
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the timescale is set by the slow relaxation time (Eq. 
[2|), leading to merger times orders of magnitude longer 
than tn- In contrast, when the MP number density is 
high enough, or the loss-cone is small enough (lower 
BMBH mass), the loss-cone is full, the merger time 
is the minimal possible, and is determined by the size 
of the loss-cone and by the dynamical time (Eq. [J). 
These conditions hold for gas-rich mergers across almost 
the entire mass range (M12 < fewxlO^MQ), and are 
also true for lower mass BMBHs in mergers of ellipti- 
cal galaxies(Mi2 < fewxlO^Afo). For higher BMBH 
masses there are not enough MPs to refill the loss-cone. 
However, the merger still evolves faster by a factor of 
~ /i2 than it would with stellar relaxation alone (table 
[3]), until the BMBH separation decreases, the loss-cone is 
filled, and the scattering rate reaches its maximal value. 
This fast MP-driven evolution continues until the BMBH 
shrinks to the point where stellar relaxation alone can fill 
the loss-cone. Since the BMBH spends most of its time 
in those late stages, the overall decrease in the total dy- 
namical merger time, tdyn, is intermediate between the 
maximal possible value of ~ fj,2^ and that due to stars 
alone (Eq. ED, 1 « t^/Jt^^^,, « 

The results indicate that MPs drive rapid coales- 
cence of BMBHs in less than tn in most minor and 
major mergers. Moreover, for most BMBHs coales- 
cence occurs in less than a Gyr, which is comparable 
to the dynamic al timescale of the galactic merger itself 
(|Barnes &: Hernquist 1992.) . Our results indicate that 
massive BMBHs (M12 > 10* Mq) in gas-poor ellipticals 
may take ^ tn to coalesce. However, these estimates 
omit processes that could shorten the coalescence time 
by an additional fa ctor of a few, such as rc-cjcction of loss 
cone stars (Milosavlievic fc MerrittI ^003; Bcrczi k et al.l 
120051 ). We find that MP-induced loss cone refilling is 
effective in driving BMBHs of M12 ~ lO"^ Mq (10* Mq, 
10^ Mq) to coalescence in a Hubble time if 0.005 (0.05, 
0.5) of the total mass density in the galactic nucleus is 
in clumped gas components with a mass function such 
as observed in the GC (Figure [T]). Since the cores of 
low-mass ellipticals with M12 < 10* Mq quite possibly 
contain some clumped gas (up to 0.02 of the MPs as- 
sumed in our major merger model; tabled]). These are 
sufficient for inducing rapid coalescence even in such sys- 
tems, even if relaxation stellar clusters or other coherent 
stellar structures is too slow. 

Figure ^ shows the evolution of the BMBH separa- 
tion for M12 = 10^, 10^, 10* and 10^ Mq in major mergers 
((3 = 1), with and without MPs. The BMBH separation 
is evolved up to the point where the decay is dominated 
by GW and coalescence follows soon after (the transi- 
tion criterion hdyn ~ ogw and Eq. [8] imply that the 
evolution curves steepen sharply beyond the transition 
point). The evolution of BMBHs with MP relaxation 
exhibits a short initial stalled phase, where the initially 
large loss-cone is empty even in the presence of MPs, fol- 
lowed by a phase of rapid decay. It should be noted that 
the decay phase does not display _the a oc evolution of 
a hard BMBH, expected when C~ const. The accelera- 
tion of the infalling stars in the Galactic potential softens 
the encounter with the BMBH and substantially reduces 
the energy extraction efficiency. Figure ([3|) shows this 
efficiency strongly depends on both the distance from 
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Fig. 2. — Evolution of the BMBH separation from as to agw in 
a major merger due to 3-body scatterings of stars. The evolution 
in the major merger scenario with MP-induced relaxation (solid 
line) is compared to that with stellar relaxation (dashed line) for 
BMBH masses of 10^, lO'', 10* and 10^ Mq (from bottom up). 
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Fig. 3^ — The dependence of the m ean dimensionless extracted 
energy Ccft = C(Q,i^eg) (Eq. IB14t for Q = 1, on the point of 
origin of the deflected star, r* /a^ , for different stages of the BMBH 
evolution a/oa (indicated by the numbers adjacent to the lines), 
taking into account the acceleration in the Galactic potentiah The 
dashed horizontal line at the top is the asymptotic value of Ccff in 
the hard limit (a/as — >0). The dash-dotted line is the value of C 
for the case a/ as = 1, when the Galactic potential is neglected. The 
vertical lines indicate the BMBH's radius of dynamical influence 
r/i and the inner cutoff of the MP distribution tmp- Most of the 
stars are deflected toward the BMBH from r* > t'mp ■ 



which stars are deflected to the BMBH and the BMBH 
separation. It should be emphasized that acceleration 
by the galactic potential cannot not be neglected, since 
it substantially reduces the efficiency of any BMBH sling- 
shot mechanism, in particular tho se wh ere the potential 
gradient is steep (e.g. [Zicr 200(| l2007f) or those where 
stars arc deflec ted to the MBH from very large distances 
(e.g. iBerczik et al.ll2006( ). 

8. IMPLICATIONS OF MP-INDUCED BMBH 
COALESCENCE 

8.1. Observations of BMBHs 
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BMBHs can be observed as resolved objects in the ini- 
tial dynamical friction stage, when a > ah, and pos- 
sibly also in the second dynamical decay stage when 
Qgw <a <ah (in particular massive BMBHs, whose an 
is large. E g. [1]). They in ay be detected indirectly (see 
review bv lKomossal[200fih . or by GW emission in the 
last GW-driven decay stage, when a < acw- Efficient 
MP-driven BMBH mergers progress rapidly through the 
second dynamical decay stage. Thus, a prediction of the 
MP merger scenario is that observed BMBHs should fall 
into a bimodal distribution: those with a > au and those 
with a < flGW) where oqw ^ a/t- In contrast, dynami- 
cal scenarios that lead to stalling, such as relaxation by 
stars alone, imply that BMBHs with a <ah should be 
common. 

The few direct and indirect observations of BMBHs 
available today are consistent with the predictions of 
the MP scenario. There are two known resolved dou- 
ble active galactic M BHs, NGC6240 with a = 1.4 kpc 
Komossa et al.' '2003*), and 0402+379 with a = 7 pc 
([Rodr iguez et al . 2006), just outside its hardening sepa- 
ration a/i~3.5 pc for Mi2~1.5 x IO^Mq. X-ray shaped 
radio galaxies, double-double radio galaxies, with pairs of 
co-aligned symmetric double-lobed radio-structures, and 
AGN with semi-periodic light curves or double peaked 
emission lines, were suggested as sig natures of clos e 
(a<Cl pc), or recently merg ed BMBHs (|Komossall2006l ). 

The detection of the GW signal from coalescing 
BMBHs would constitute direct evidence of such events. 
Our calculations show that for most galaxy mergers, the 
BMBH would coalesce within tn, and so the BMBH co- 
alescence rate should follow the galaxy merger rate. In 
that case the cosmic rate of these GW events could be 
as high as f O^yr "^ ()Haehneltl [l99l iSesana et al.l 120041 : 
lEnoki et al.l[200l . 

8.2. Triple MBHs and MBH ejection 

The galaxy merger rate in dense clusters may be high 
enough (> 10~^ yr~^; FMamon 2006) for a second merger 
to occur before the first BMBH coalesces. This would 
result in the formation of an unstable triple MBH sys- 
tem, which will eje ct one of the MBHs at high velocity 
(|Saslawet al.l[l97i) . This scenario was suggested as a 
possible solution for the stalling problem, as the third 
component may drive the BMBH to hi gh eccentricities 
and to much more r apid coalescence (Blaes et al.ll2Q02l : 
llwasawa et al.|[2006t [Hoffman fc Loebll2007D . 

Because MPs accelerate BMBH coalescence in most 
mergers, our results indicate that triple MBH systems 
and high-velocity MBHs ejected by the slingshot mech- 
anism should be rare^ (Fig- [J)- Because of the rapid 
BMBH decay, in those cases where a triple MBH is 
formed, it is expected that it will be hierarchical. This 
would t ypically lead to fast coalescence of the inner 
BMBH jMakino fc Ebisuzakilll994l:llwasawa et al.ll2006[ ). 
followed by the MP-driven decay and coalescence of the 
newly formed central MBH with the outer MBH. 

8.3. Mass deficits 

Recent observations of an apparently host-less quasar 
[Magain ct al. 2005) were interpreted as an ejec ted MBH 
fflaehnelt et al.l i200a IHoffman fc Loebl [200^. |2003), but see 
[Merritt''et'^Lr i r2006l) for an opposing view. 
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Fig. 4. — The influence radius in the post-merger galaxy at 
the stalling time, r^, as function of Q, for initial Dehnen density 
profiles with 7 = 0.5, 1.0, 1.5 (top to bottom), as derived bv lMerrit^ 
(l2006i ) in A'^-body simulat ions (crosses), and by the approximate 
analytical expression Eq. I IA5I I (lines). 



The large number of stars ejected from the sys- 
tem during the BMBH coalescence could change the 
stellar distribution of the BMBH environment. It 
has been suggested that the mass deficit observed in 
some b right elliptical galaxies is the result of such 
events (Milosavlievic et al. '2002b'; 'Ravi ndranath et al.l 
|2002; Graham 2004; Ferrarese et al. 200^). The total 
mass of ejected stars in the dynamical decay phase de- 
pends only on the initial and final BMBH separations, 

Mej(0 = / dt' f dE:F{E, t') ~ JMi2 In ^ , 

(12) 

where T is flux of stars supplied to the loss cone,_and J 
is a numerical factor approximately equal to 1/2( 7 (Q96; 
iMilosavlievic fc Merrittil20"03l : ISesana et al.ll200l . 

Previous studies of the ma ss deficit 

(Milosavljevic et alj l2002at iMerrittI 120061 ) took into 
account only the stars evacuated from the core before 
the BMBH stalled at a ~ a?i because of inefficient stellar 
relaxation. We note that that there are between 2-7 
further e- foldings between at and acw (Fig- [2])- As a 
result, when the BMBH merger is driven all the way to 
OGw, the mass deficit will grow substantia lly (this was 
recently also noted bv iMerritt et al.l[2007[ ). When the 
merger is driven by MPs, the mass deficit will be on 
the ^ (l-2)rMP scale, where most of the scattered stars 
originate. We calculated the total mass of stars that 
originated from such distances and that were ejected 
from the core during coalescence. We found that these 
amount to approximately 30-40% of the total stellar 
mass in these regions. 

We note that generally, the magnitude and spatial scale 
of the mass deficit could, in principle, discriminate be- 
tween different proposed solutions for the stalling prob- 
lem. In mergers driven by non-axis ymmetric potentials , 
most stars originate from large radii lBerczik et all ([20061 ) 
where the enclosed number of stars is very large. The 
fractional, spatially-averaged mass deficit will therefore 
be very small, and harder to detect. In contrast, sce- 
narios that assume very steep cusps ([Zieil l2006l l2007l) 
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lead to substanti al central mass dep l etion. Even in ga s 
induced mergers (jEscala et all 120051 : iDotti et al.ll2007f ). 
where stars play a minor role, there may be an indirect 
mass deficit effect caused by the inhibition of star forma- 
tion due to heating of the gas by the inspiraling BMBH. 

9. DISCUSSION AND SUMMARY 

We have shown that it is very likely that MPs play a 
dominant role in the aftermath of galactic mergers. They 
shorten the relaxation timescale in the galactic nuclei by 
orders of magnitudes relative to 2-body stellar relaxation 
alone, and drive the newly formed BMBH to rapid coa- 
lescence by the slingshot effect. The MP mechanism re- 
quires only the existence of large enough inhomogencities 
in the galactic mass distribution. Since these occur natu- 
rally over a wide range of conditions, the MP mechanism 
is robust and probably accelerates most BMBHs merg- 
ers. The one possible exception may be mergers of two 
gas-poor elliptical galaxies, where GMCs are less com- 
mon. However, simulations indicate that stellar clusters 
can play the role of MPs and drive an efhcient merger 
even in most of these cases. 

We conclude that most BMBHs are expected to coa- 
lesce within tn, even in cases where previous theoreti- 
cal modeling, which did not consider accelerated relax- 
ation by MPs, predicts that the merging BMBH stalls. 
This conclusion is based on conservative assumptions. 
We considered only circular BMBHs, whereas eccentric 
BMBHs coalesce even faster in the GW emission dom- 
inated phase, and we neglected the possible concurrent 
effects of any of the other orbital decay mechanisms pro- 
posed in the literature. It thus appears likely the BMBH 
coalescence is in fact achieved on timescales ^ tH , which 
implies that BMBH coalescence GW events occur at the 



cosmic rate of galactic mergers. Specifically, we predict 
that low-mass BMBHs, which are prospective LISA tar- 
gets, should coalesce within a merger dynamical time, 
lO^-lO^ yr. 

Efhcient MP-driven BMBH coalescence have addi- 
tional implications, which were discussed here briefly. 
Fast BMBH mergers decrease the probability of nuclei 
containing triple MBHs, and hence of ejected MBHs, 
since in most cases, the BMBH coalescence time is 
shorter than the mean time between galactic collisions. 
During the final stage of the merger, when the BMBH 
separation shrinks from the hardening radius to the fi- 
nal GW radius, a large number of stars will be ejected 
from the nuclei. We find that this additional ejection 
stage could appreciably increase the mass deficit of the 
newly formed nucleus, beyond what is predicte d takin, 
into a ccount only the earlier stages of the merger (|Merrit 
l2006h . 

In summary, we have shown that the plausible exis- 
tence of MPs in galactic nuclei shortens the relaxation 
time by orders of magnitude. In particular, MPs ac- 
celerate the dynamical decay of BMBHs by efficiently 
supplying stars for the slingshot mechanism. This pre- 
vents stalling (the "last parsec problem" ) and allows the 
final coalescence of the BMBH by GW emission within 
a Hubble time. 
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APPENDIX 
A. THE STALLING RADIUS 

This appendix presents a simple analytic approximation for the stalling separation, a^, as function o f the pre -merger 
galactic density profile and the BMBH mass ratio Q, which is based on the A^-body simulations of iMerritH ([2006). 
Typically, ^ (Eq. [1]), up to a factor of a few. Assuming the ansatz — s- in the evaluation of the BMBH 
coalescence time can lead to inaccuracies of up to a factor of a few, in particular for Q ^ \. 

Merritt ( 20 0^ model ed typical galactic cores in large A^-body simulations of BMBH coalescence by Dehnen config- 
urations (|Dehnenlll99"3l) . 

M 1 

[4V(3-7)]d3 (r/d)'^(l-Hr/d)4-7 ' ^ > 

where M is the total stellar mass, d a scale length and —7 the logarithmic slope at r ^ d. A central MBH of mass 
Mi/M— 0.01 was added to the initial density distribution. We assume here that the results derived for this particular 
class of models also apply, at least approximately, to other initial density distributions and MBH-to-stellar cluster 
mass ratio s. 

iMerrittI ()2006l ) finds that the stalling radius can be described to a good approximation, independently of 7, by 

as = 0.2Q/(1 + Q)^r'^iMi2) = 0.8 K(Mi2)/r,,(Afi2)] an , (A2) 

where r^(Mi2) is the radius of influence of the BMBH at the stalling time tg, after the scouring effect of the binary 
formation, which is estimated as follows. The ejected mass at tg can be approximated analytically. 

' . (A3) 



M- 
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The post-merger radius of influence rj^ can be estimated to better than 3% typically (Fig. 0]), by assuming that the 
post- density profile resembles the original profile, except for the removal of AM from the center further out, so that 



M(< r;j = 2Af{2 = 2A/i2 + AM = Mi(l + Q){2 + O.Tg"-^) . 
The enclosed stellar mass in the initial Dehnen distribution is M{< r) — Ad[r/{r -f d)Y^'^ ^ and so 



r;(Mi2)/r;,(Mi2) - (2A/12/M) 



a/(7-3) 



- 1 



(2M{2/M)'/'^^^) - 1 



(A4) 



(A5) 
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The correction is thus a function of the inner cusp slope only. The stalling separation (Eq. IA2p is a rising function of 
Q. For Q^O, as-^Q.Sah irrespective of M12/M or 7. For Q ~* 1, M12/M = 0.01 and 7 = 2, as'^2.2ah- 

B. BMBH ENERGY EXTRACTION BY INTERACTIONS WITH STARS 

This appendix details how the mean BMBH energy that is extracted by an encounter with a star from the galactic 
core is estimated using results of isolated 3-body scattering experiments, which are available in the literature. 

B. 1 Adaptation of results from scattering experiments 
The rate at which the BMBH changes its binding energy £'12 — GMiM2/2a due to interaction of stars is 

{AE{b;Q,0)2nb—db, (Bl) 



dt Jo db 

where b is the impact parameter, ^ = 1^/^12 is the hardness parameter {v is the velocity of the incoming star at infinity 
relative to the BMBH center of mass, ignoring the potential of the galaxy) and V12 = ^GM^/a is the BMBH's circular 
velocity in the reduced mass system. The quantity {AE(b; Q,C)) is the mean energy extracted by the star from the 
BMBH orbit averaged over the BMBH orbital parameters, and 27rb(dr/d6)d6 is the rate at which stars are deflected 
into orbits with impact parameter in the range {b, b + db). In a Keplerian system the impact parameter b = x6o, with 
bQ = 2GMi2a/v'^ = 2a'^/£,^ (Q96) is related to the periapse distance rp = ya by b'^=rp (l + 2GMi2/rpV^) , which can be 
written as 

x'^eyV2+y, (^^i+2ex^-i)/ e- m 

The extracted energy (A£') is a function of Tp, the mass ratio Q and the hardness parameter £. It was derived 
numerically by Monte Carlo simulations of 3-body scattering ( iHilli [19831 : iQuinlanI 119961 : ISesana et al.ll200^ . which 
sho w that it i s large and fairly constant for rp/a < 0.5 — 2, and then falls rapidly to zero for rp/a > 0.5 — 2 (see 
also IZieij|2007l for an extended discussion of the rp-dependence of the extracted ene rgy). Unfort unately, the behavior 
of {AE(b]Q,(,)) as function of its parameters has only been partially documented. iHillsl ()1983l) stud ied t he 6 and Q 
dependence only in the ^ ^ limit, for specific values of the BMBII orbital eccentricities, ISesana et al.l (j2006l ) show 
plots only for ^ ^ 0, while Q96 explored the full range of values for ^ and Q and averaged over the eccentricity, but 
integrated over the b dependence. For that reason, it is not possible to use these results to evaluate Eq. (jBl[) explicitly, 
and it is necessary to resort to an approximate formulation. Following the trends seen in the ^ — > simulation results, 
we adopt here a step function approximation, which is based on the assumption the (AE) is roughly constant between 
6 = and an effective maximal impact parameter 61, and write 



dEi2 
"df 

where the ^-averaged extracted energy is defined as 



AEiT{bi), (B3) 



AEi = {AE{Q, 0) = ^'^ 2iib {AE{b- Q, 0) db j nbj , (B4) 
and where the total rate of stars with impact parameter b<bi is 

/■^i dr 

r(6i) = / 2TTb—db. (B5) 
Jo d& 

It should be noted that the step function approximation implies that the periapse-averaged energy extracted per star 
should not depend strongly on the mode of loss-cone replenishment, whether it is in the empty loss cone regime, where 
stars diffuse into the loss-cone from its boundary (r^ ~ a), or whether it is in the full loss cone re gime, where the 
orbits span the entire range 0<rp<a. This behavior is indeed seen indirectly in the simulations of lMerritt fc SzelH 
()2006l Fig. 2), where the mass ejection rate remains nearly constant as the system transits from the full to the empty 
loss-cone regime. 

Q96 does not quote the simulation results in terms of (AE) directly, but rather in terms of a related dimensionless 
quantity Hi, which expresses the rate at which energy is extracted from the BMBH by an ambient background of 
stars with mass density p and velocity v at infinity, d{l/a)/dt — HiGp/v. IHillsl (fT983i ) and Q96 express the extracted 
energy in a dimensionless form, 

C= M12AE/ 2AUEi2^aAe/Gp , (B6) 

where Ae = AE/Mi, is the specific energy extracted by the star and p = A/1M2/M12 is the reduced mass (C = 1 
corre spo nds to the case where the specific energy carried by the star equals twice that of the BMBH). In terms of C, 
Eq. (|B4p can be written as 

Ci= r 2nb {C{b- Q, 0) db / TTbl (B7) 
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The quantity Hi is related to the orbitally averaged (C) in terms of the dimensionless impact parameter x by (Q96, 
Eq. 11) 

Hi{Q,^) =8tt x{C{x;Q,£,))dx^ 4:TT J 2TTb{C)db^ TTbl- (B8) 



A comparison of Eqs (|B7p and (jBSp shows that Ci is related to the values of Hi quoted by Q96 through the relation 

CiiQ,O^Hi{Q,0{bo/bi)y47r. (B9) 
Let rpi = yia be the periapse that corresponds to the effective impact parameter 61. It then follows that 

CiiQ, = HiiQ, 0/ [47ryi(l + yie/2)] . (BIO) 
Q96 suggests an analytic fit to the behavior of Hi as function of Q and ^ (Q96, Eq. 16 and table 1), 

HiiQ,0= ^ ""'^^^ , (Bll) 

V'i + KM(g)f 

where here we use a slightly different notation from Q96, wo = WQgG/Vi2, so that the values of wo{Q) are the numeric 
values in the 3rd column of table 1 in Q96, or given by the analytic fit formula (Q96, Eq. 17) 

wo{Q) ^ 0.85 V^2/A/i2 = 0.85 + Q) ■ (B12) 

The tabulated values of Hq[Q) (Q96, table 1, 2nd column) are nearly independent of Q, with Ha = 21.lt_W over the 
range Q = 1/256 to Q = l (the lower values corresponding to larger Q). 

For a distribution of velocities characterized by a ID velocity dispersion cr, the effective H is given in terms of 
C^<^/Vi2 by (Q96, Eq. 18) 

H{Q,0=Hi(Q,Vic)[VV^ + ^og[l + a[awo{Q)f]] , (B13) 
where a = 1.16 and /3 = 2.40. This can then be approximately related to the mean ejection energy as in Eq. (jBlOp by 

C(g, C) = H{Q, 01 [4^2/1(1 + yiCV2)] . (B14) 

The nature of the approximation here is that the translation between b and is done for a representative velocity, 
assuming Keplerian motion. For the purpose of numeric calculations we assume yi = l. 

The relevant scale for MPs is rMP^2r/i, which encloses >4M, in stars. On that scale the potential is dominated by 
the stars. For a stellar density distribution far from the MBH, the velocity dispersion is (T^(r) ~ G'Mi,(< r)/r~ 
const. Here we represent the typical initial stellar velocity by the circular velocity v = y/2a. This velocity needs to be 
corrected for the fact that the star is accelerated by the galactic potential as it falls toward the MBH. In the fictitious 
3-body syst em ( the BMBH and the star), its effective hardness parameter , (lcff{r)>C,, depends on the star's point of 
origin (see i]B.2[) . Thus, the BMBH total decay rate is given by integrating over the contribution of stars originating 
from all radii, with the loss-cone size expressed in terms of the periapse, 

- -2—— / C[Q,C,cs[r;a)\ dr. (B15) 



d< M12 J ^-^'•'"'^ ' " dr 

Note that the "hard limit", d(l/a)/dt = const, is recovered when the cross-section for interacting with the BMBH is 
dominated by the gravitational cross-section term (C 0, C = const, Eq. IB14p . and Fcxa (full loss-cone regime, Eq. 

El). 

B.2 Dependence on the galactic potential 

The extracted energy (Eq. IB14|) depends on the hardness of the encounter C- In the soft encounter limit (C ^ 0), the 
interaction with the BMBH is strong and independent of C, (Eqs. IB10| IbTT|) . The hardness of the encounter depends 
on the star's point of origin. The farther away it starts from the center, the more it will be accelerated by the galactic 
potential, and the faster it will be when it approaches the BMBH. Since MPs typically deflect stars into the loss-cone 
from large distances, the effect of the galactic potential in making the encounters harder and less efficient cannot be 
neglected (Figure [3]). 

The 3-body scattering experiments of Quinlan (1996) took into account only the potential of the BMBH. To relate 
the energy extraction to that of a star falling in the combined potential of the BMBH and the galaxy, it is necessary 
to calculate the corresponding effective initial velocity the star should have in the fictitious 3-body system containing 
only the BMBH and the star. 

The gravitational potential in a spherical stellar system with mass density p{r) surrounding a central BMBH of mass 
M12 is 



r) = -AttG 



p{x)x dx + p{x)xdx 







(B16) 
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For a power-law mass density profile between ri <C ^2 , 



Pir) = Pa[ — 
Jo 



the enclosed stellar mass is (a 7^ 3, r^ri) 



M(< r) 



47r 
3 — a 



3-a 



The potential at ri < r < r2 is 



0(r) = ii_W _ , p(x)xdx 



GMi2m{r) 



47rG'/9o?'o 



Q / ln(r2/r), 



a = 2 



(B17) 



(B18) 



(B19) 



where m(r) = 1 + Af(<r)/A/i2 is the total mass up to radius r relative to the BMBH mass. 

Suppose a star starts falling toward the BMBH with velocity v from an initial radius r down to a radius r\ close 
enough to the BMBH so that M(< ri) ^ 0. The specific orbital energy of the star, e, is conserved, 



(B20) 



Taking the velocity of this star at as the velocity it would have if it began falling from the same distance in the 
fictitious 3-body system, we can find what the effective velocity, Ueff, it should have at r far away from the BMBH 



where 



The effective velocity it then 



4>\r) = 



GMi2 



Is = 1-2 + 2W{n) ~ 4>'{r) + </)(r) - 0(ri)] = i;^ + 2GAfi2 



1 1 niir) mirA A'kG 

\ \-L^ ^—^ 

Ti r r Ti M- 



•J 12 



p{x)xdx 



Since rn(ri) — > 1, 



2G 



M{< r) 



+ An p{x)dx 



For the a = 2 power law assumed here (r^ri). 



(B21) 
(B22) 

(B23) 
(B24) 

(B25) 



The effective hardness parameter is then Ccff (^) = Wcff/Vi2, which is used to evaluate the mean extracted energy (Eq. 

[m. 
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